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Abstract 

We study many-body correlations in the ground states of a general quantum system of bosons or 
fermions by including an additional Jastrow function in our recently proposed variational coupled- 
cluster method. Our approach combines the advantages of state-dependent correlations in the 
coupled-cluster theory and of strong, short-ranged correlations of the Jastrow function. We apply 
a generalized linked-cluster expansion for the Jastrow wavefunction and provide detailed analy- 
sis for practical evaluation of Hamiltonian expectation value as an energy functional of the Jas- 
trow function and the bare density-distribution functions introduced and calculated in our earlier 
publications; a simple, first-order energy functional is derived and detailed formulas for higher- 
order contributions are provided. Our energy functional does not suffer the divergence as in most 
coupled-cluster calculations when applying to Hamiltonians with hardcore potentials. We also 
discuss relations between our energy functional and the energy functionals from other theories. 

PACS numbers: 03.65. Ca,31.15.Dv,31.15.Pf,71.10.-w,71.45.-d 
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I. INTRODUCTION 



Most microscopic quantum many-body theories developed over the last five decades can 
perhaps be broadly divided into two categories, one in real space and the other in momen- 
tum (configurational) space. While a real-space many-body theory usually focuses on the 
interaction potential part of Hamiltonian and evaluates the Hamiltonian expectation value 
in a first quantization form, a momentum-space theory often starts from the kinetic part of 
Hamiltonian, is a basis of or closely related to many-body perturbation theories, and mostly 



deals with Hamiltonian in a second quantization form 



A typical real-space approach 



to the ground state of a quantum many-body system is prov.ded by Jastrow^wavefunetion 
which is constructed by a state- independent two-body correlation function [3]. Systematic 
techniques based on Jastrow wavefunction including extension to inhomogeneous boson or 
fermion systems is now generally referred to as the method of correlated basis functions 
(CBF) . The CBF method has proved to be efficient in dealing with strong, short- 

ranged correlations typified by those in quantum helium liquids jj]. On the other hand, 
momentum-space many-body theories are often easier to apply and, due to inclusion of the 
state-dependent correlations, is capable of producing accurate results for a wide range of 

n " 



quantum systems, such as boson gas [7], quantum antiferromagnets with Neel order 
finite nuclei Isl, and electron systems such as electron gas [10], atoms and molecules A 
typical momentum theory is the coupled-cluster method (CCM) in which wavefunction are 



explicitly constructed by state-dependent operators [12|, ll3||. State-of-the-art calculations of 
the CCM with high accuracy have often been carried out in quantum chemistry [l^, and 



recently in quantum spin lattices with Neel order [15|. Systematic resummations of dia- 
grams in perturbation theory for boson systems have revealed interesting relations between 
the two approaches, for example, the hyper-netted chain approximation in the CBF method 
in fact contains a consistent resummation of both infinite ring and infinite ladder diagrams 



of momentum-space approach 16|, \V£\. It appears that real-space and momentum-space 



approaches complement each other and unification of these two approaches may provide a 
quantitative description applicable to wider range of quantum many-body systems, including 

n 

in particular the strongly-correlated fermion systems [18]. 

We recently extended the CCM to a variational formalism in which bra and ket states are 



now hermitian to one another 



3, 



2Q|, 



211 ] ■ contrast to the traditional CCM where they are 



2 



not 



22|. We introduced the hermitian-conjugate pair of important bare density distribution 



functions for practical and systematical calculations; the traditional CCM was shown to 
correspond to a simple linear approximation in one set of distribution functions of our 
variational coupled-cluster method (VCCM). The well-known momentum approaches such 
as Bogoliubov theory of boson gas, Anderson's spin- wave theory (SWT), and BCS theory 
of superconductivity 23j, can all be explicitly shown as special low-order approximations 
in both the ground- and excited-state wavefunctions of the VCCM. We have demonstrated 
by a detailed application to quantum antiferromagnets with Neel order. Approximations 
beyond SWT improved results for the ground-state properties and new excitation states 



have also been obtained 



2l| . Furthermore, our calculations for the bare density distribution 



functions can be carried out by diagrammatical techniques similar to those employed by 
the CBF methods {2^. Hence, a bridge between coupled-cluster theory and Jastrow theory 
is built. We therefore believe it is a natural next step to combine the two methods for 
a unified description. This is our main purpose in this article. Krotscheck, Kiimmel and 
Zabolitzky (KKZ) made the first attempt in 1980 for fermion systems ISj. They employed 
the traditional CCM with only the ket state specified and the Jastrow function is fixed before 
hand; Hamiltonian eigenequation was used to obtain the ground-state energy and equations 
for the ket-state coefficients. Here we employ explicit ket and bra states of the VCCM and 
calculate distribution functions in terms of these ket- and bra-state coefficients; the energy 
functional is derived in terms of the distribution functions and the Jastrow function. 

We focus here on formal development of our approach and organize this article as follows. 
In Sec. II we introduce our correlated wavefunctions and their distribution functions. The 
generalized linked- cluster expansion technique is employed for calculations of generating 
functional in terms of the bare distribution functions. In Sec. HI we evaluate Hamiltonian 
expectation value using Jackson-Feenberg transformation; a first-order energy functional is 
derived and formulas for higher-order contributions are provided. We conclude in Sec. VI 
with a summary and discussion on the relations between our approach and other many-body 
theories. 
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II. DISTRIBUTION FUNCTIONS IN VARIATIONAL JASTROW COUPLED- 
CLUSTER THEORY 



As proposed earlier [20|, we consider a general variational wavefunction by adding a 
Jastrow correlation operator on top of the Coester state \'^'^) 



\^) = e^°/^\^'), (1) 



and we use the following notation 



(^1 = (^-|e^°/^ (2) 



for the bra state. In Eqs. (1) and (2), the hermitian operator is the Jastrow correlation 
operator, and is given in terms of field operators as j^, 

5'° = i /" dxidx2ip^{xi)ip^{x2)u{xi,X2)ip{x2)ip{xi), (3) 



2 

where field operators and ip{x) obey the usual boson or fermion commutation relations, 
u{xi,X2) is a local, symmetric function u{xi,X2) = u{x2,Xi), and x are particle coordinates 
including spin degrees of freedom. We also require that u{xi,X2) is bound and short-ranged 
in real space. We do not include in any single-body operator as it can be easily absorbed 
in the Coester states to be defined next. The boson or fermion symmetry is contained in 
the Coester states by definition. More specifically, the Coester ket state is given by, using 



the convenient notation invented by Arponen and Bishop 24l |. 



|^<^) = e^|<l>), (4) 

where S is constructed by the so-called configurational creation operators Cj which are 
defined with respect to the model state |$) with the nominal index / labeling multi-particle 
excitation states from the model state |$), 

^ = Ewl, (5) 

with Fj often referred to as the correlation coefficients. The Coester bra-state (\E''^| is given 
by the hermitian conjugate of the ket state, 

(^^| = ($|e^ S = Y.FiCi (6) 

I 



in our VCCM, where Cj are the corresponding configurational destruction operators and Fj 
are the independent, hermitian conjugates of F/. In the traditional CCM, however, the bra 
state is parametrized differently from the ket state and is written as [22] , 



(^'=1 = ($|^'e-^, ^' = 1 + ^ F'jCi. (7) 

where S is defined as in the ket state of Eq. (5), and F[ is the bra state coefficients which, 
in general are not hermitian conjugate of Fj. As discussed in Appendix A, the CCM states 
violate the condition for application of the generalized linked-cluster expansion due to the 
linear construction of the bra state. Furthermore, the evaluation of the kinetic energy dis- 
cussed in Appendix B will contain a three-body terms because the nonhermitian relation 
between the ket and bra states in the CCM. We therefore will not discuss the CCM further 
and focus only on the VCCM basis of Eqs. (4-6). Clearly, the natural variational parameters 
are {F,F,u) in the VCCM basis, where we have used the notations F = {Fj}, F = {Fj} 
and u = u{xi,X2)- In principle, if these Coester states are exact (namely, all configurations 
are included in the summations over all /-indices), parameter u is redundant. However, as 
we always need to make a finite truncation approximation in summations over /-indices in 
any practical application and it is well-known that the Coester states in a finite trunca- 
tion approximation are not efficient in dealing with the strong, short-ranged correlations, 
the two-body Jastrow function u{xi,X2) is a useful, important variational parameter in a 
real application. We hence always assume the summations in Eqs. (5) and (6) are within 
the subset of a truncation approximation. One of such truncations is the so-called SUBm 
approximation in which we retain up to m-body creation operators only. 

Our basic strategy for calculations is to evaluate the generating functional W of Eqs. (1) 
and (2), 

= In(^'l^) = M/'^ + M/", (8) 
where is the generating function for the pure Coester states without the Jastrow operator, 

W = ln(^"|^"), (9) 

and is the remainder containing the u function 

W- = ln{e'y, (e^«)c=iS^!>l. (lo) 
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The general strategy for calculating = W^i^F, F) was discussed in our earlier papers and 
detailed calculations were demonstrated in the spin-lattice application [l^, Q. Briefly, we 
first introduce the hermitian-conjugate pair of the bare density distribution functions 

9i = (c\r = ^^ ^?/ = (c^/r = ^, (11) 

where the expectation values (■ ■ are calculated using the Coester states of Eqs. (4-6) 
as was given in Eq. (10); any physical quantity is expressed in terms of these distribution 
functions; these distribution functions are then calculated either by algebraic technique |l^ 



or by diagrammatic technique 20|] . In the algebraic approach we derive the following general 



and self-consistency equations for the calculations, 

ri = GjiF,g^), g^j = Gj{Frf) (12) 

where Gi{F,g'^) with notation g'^ = {gj} etc. as before is a function containing up to linear 
terms in g'^ and finite-order terms in F, and similarly for Gi{F,g'^). In the diagrammatic 
approach, many techniques employed by the CBF method are also applicable here. For 
example, the coefficients Fj can be replaced in every diagram of expansion by the bare 
distribution functions gi after resummations of infinite ring diagrams. In general, the energy 
expectation value can be expressed as a functional polynomial in F and g'^, E'^ = E'^{F,g'^); 
and the variational equations can then be derived by the functional derivatives [zS]. 

In terms of the distribution functions as discussed above, we now consider evaluation of 
{e^°)'^ in Eq. (10). After cluster expansion for the Jastrow part of the wavef unctions, we 
apply the generalized linked-cluster theorem discussed in Appendix A in a similar fashion 



261], and we obtain 



as the traditional Jastrow theory 

= ln{e'"r =T.^ fdxi- --dxnlZnPlU (13) 

n=2 

where Zn{Y) are the n-body Yvon-Mayer functions of the bound function Y12, 

= e"(^i'^2) - 1; (14) 

is the ra-body density distribution functions of the Coester states; and the notations [■ ■ -Jl 
denote that only the linked components of the products are included. In deriving Eq. (13) 
in Appendix A, we have assumed that the Coester states satisfy the cluster-decomposition 
properties (i.e., higher-order are given by sum of products of lower-order ones plus a 



nondecomposable core). While this is true in a SUB2 approximation in our spin lattice 



app lication and the SUB2 approximation for the Bose gas and BCS superconductors 



271 ]. we can not provide a general proof without the details of truncation approximations 
employed in the Coester states in a real application. Hence by using Eq. (13), we have made 
a further restriction in the truncation approximations of the Coester states. In the following 
discussion this further approximation is understood and we will examine this property in 
real applications beyond the SUB2 approximation. The following analysis for calculating 
in terms of Fj and Qj clearly show such examination posing no major difficulty. 
As all the Yvon-Mayer diagrams in Z2 and Z3 are linked, we simply have 

Ml = l^i2P2, P2(^i,^2) = {i^\x,)i;\x2mx2Mx,)y (15) 
for 2-body cluster contribution and 

[ZsPtU = {Y12Y2S + 1^23^31 + 131^2 + Yi2Y2sYsi)pl (16) 

for the 3-body cluster contributions with the 3-body distribution function 

= {4j^{x,)^^{x2)tPHx3)4j{x3Mx2)4j{x,)y (17) 

in the Coester states. From the 4-body function Z^iY) and onward, however, there are 
unlinked Yvon-Mayer diagrams which are to be included in the product [Znp'^^L only after 
multiplying with the terms of to form the linked components, and details of which will 
depend on applications with a truncation approximation in the Coester states. There is 
clearly a trade-off between the order of linked-cluster expansion and the order of truncation 
approximation of the Coester states. We hope to get experience in real applications in future. 
The additional state-dependent correlations of the Coester states are contained firstly in 
of Eq. (9) and secondly in the density distribution functions p^. Importantly, these density 
distribution functions can be calculated in terms of F and g'^, and by using the linearity 
theorem of the VCCM j21|, we can show that all p^ functions contain only up to linear terms 
in g'^ and finite-order terms in F. As a demonstration, we consider the two-body function 

Pi. 

Pi(xux.) = (A,Y = ■jJ^^mA.e'm = ^^(¥\e'Am, (18) 

where A2 = ip'^ {xi)ip'^ {x2)4'{x2)i'{xi) and A2 is 

A2 = e-^A2e^ = e-^'ijj\xi)ij^{x2)tlj{x2)ij{xi)e^ = i/j^ {xi)tlj^ {x2)tlj{x2)ij{xi), (19) 



with (x) = e ^ilj'^[x)e^ etc. Using the nested commutation series, we obtain 



^\x) = i:\x) + [ij\x), S] + -[[ij\x), S],S] + 



^\x) + [^\x),S], 



(20) 



where the series terminates at first order as S is constructed by creation operator Ct only 

nn 

|12l. Il9l|. Evaluation of v42|$) in general leaves only a constant and creation operators acting 
on 1$), namely, 



X2,o{F; xi, X2) + X2j{F; Xi, X2)Cj 
I 



(21) 



X2fi and X2J are a two-body function containing up to fourth-order terms in F. Therefore, 
the two-body density function in the Coester states is, using the definition of Eqs. (11) 

P2{xi, X2) = X2fi{F; xi, X2) + X2,i{F; Xi, Xs)^/. (22) 

I 

In similar fashion, we derive 

pl = Xn,o + T.Xn,i~g'j (23) 
I 

for the n-body density distribution function in the Coester states, where Xnfl and Xn,i are 
the ra-body functions containing up to (2r2)th-order terms in F. Therefore, the linked-cluster 
contributions of Eq. (13) is written as 



W=- I dXidX2Z2 



X2fl + ^ -^2,/^/ 
I 



H — / dxidx2dx3Zs 
6 



^3,0 + Y 



+ 



(24) 



where the remainders are the 4-cluster and higher-order contributions, and their calculations 
will depend on the details of the truncation schemes employed in the Coester states but the 
general property of finite-order in F and linear in g'^ remains. 

Before we consider the density-distribution functions of the Jastrow-Coester states of 
Eqs. (1) and (2), we need to define biased distribution functions as 

gi = ^^me^'cl\^'), gi = ^^(mCie^'\^'). (25) 

They are so called because they are not defined usually as (C|) and clearly gj 7^ (C|) 
due to the fact that C] and S'°/2 do not commute in general. Similarly, gj 7^ (C/) by our 
definition. These biased distribution functions can be calculated by the functional derivative 
of the generating functional of Eq. (8) as 



9i 



dW 



~9'i + 



(26) 



where we have used the definition of Eq. (11). The functional derivative in Eq. (26) can be 
calculated by Eq. (24), 



-gP^+X2,i'{gi+r-gigp) 



+ ■■■, (27) 



where we have used the fact that 

with gj^j, = {ClC\,y = {C\^i,y. Clearly gj^P are also bare distribution functions. 

Using the fact that commutes with density operator ijj^{x)ijj{x), the single-particle 

density function, pi{x) — ('0^(x)'0(a;)) can then be calculated as 

pi{x) = -^J^{-^^\e^°ilj\x)ilj{x)e^\^) = -jJ^{-^^\e^'e^ij^{x)ilj{x)\^) (29) 

where the evaluation of -j/^^-i/^l^) = (Xi^o + J2i Xi,iCj)\^) is similar to that in Eqs. (19-23), 
hence we have, using the definition for the biased distribution function of Eq. (25) and then 
using Eq. (26), 

p,{x)^p'i(x)+j:X^j{F;x)^, (30) 

where pi {x) is the one-body density-distribution function of the Coester states as in general 
given by Eq. (23) and the functional derivative dW^/dFj is given by Eq. (27). 

The two-body density distribution function of states of Eq. (1) and (2) can be calculated 
in similar fashion as pi shown above. We take a more efficient calculation by the functional 

derivative as 

dW dW 
= = 25-^ = 2^-^. (31) 

The functional derivatives dW^/ du only involve the Yvon-Mayor functions and it is easy 
to derive from Eq. (13) 

P2{x,,X2) = e"(^--^)p^(xi, ^2) + ^ / dx\dx'2dx'f-^^^^^pl{x',, 4, 4) + • • • , (32) 

where are the n-body density-distribution function in the Coester states as given in 
general by Eq. (23). We can see immediately from Eq. (32) that the short-ranged correlation 
function w(xi, X2) will play an important role for applications to strongly-correlated systems, 
where the pure Coester states are known to be inefficient. 

9 



In summary, we have calculated the density-distribution functions as functionals of 
(F, g'^, u) where are the bare density distribution functions discussed in our earlier VCCM 
papers these calculations are all straightforward up to and including the third-order 

cluster contributions as these functionals are all polynomials of F and g'^, and u enters into 
these functionals through Yvon-Mayor functions. The fourth- and higher-order cluster con- 
tributions will dependent on the details of applications with the truncation approximations 
employed in the Coester states. The first few terms of one-body and two-body functions 
are given by Eqs. (30) and (32) respectively. We also want to point out that evaluation 
of (C|) or {Ci) (or the off-diagonal density functions such as {■ip'^{x)ip{x'))) in general are 
highly nontrivial due to the fact that these operators do not commute with the Jastrow fac- 
tor exp(S'°/2). However, these calculations are not needed in the Hamiltonian expectation 
value to be discussed in the next section; if necessary these calculations can be carried out 
in an approximation after we have obtained solutions to the variational equations 28|. We 
will discuss these calculations in future. 



III. EVALUATION OF HAMILTONIAN EXPECTATION VALUE 

In evaluating a general Hamiltonian expectation value, we first notice that the kinetic 
part of Hamiltonian in general does not commute with the Jastrow operator in our 
states of Eqs. (1) and (2). In real space, however, the kinetic operator contains only second- 
order derivatives in particle coordinates. We want to take this advantage by expressing our 
states in real space for calculations. As shown in Appendix A, the wavefunctions of the 
Jastrow- Coester states of Eqs. (1) and (2) in real space are given by a product 

^ = ^'^^^^ ^ = (33) 

where = e^^^ = exp[Y,i^j u{xi, Xj) /2] is the familiar Jastrow wavefunction and 
'^'^{xi, ■ ■ ■ ,xn) and \E'^(xi, ■ ■ ■ , xat) are real-space wavefunctions of the Coester ket- and 
bra-states respectively. In general, we do not need to know the explicit functional form of 
and as our calculations involving them are always carried out in a second quantiza- 
tion form as we show below. It is interesting nevertheless to know that in a low-order SUB2 
approximation, many-body function is known explicitly as a partial-wave function, the 



so-called independent pair functions for boson gas 



29| or BCS superconductors [30]. The 
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Coester wavefunctions and obey proper symmetry, namely they are antisymmetric 



for fermions and symmetric for bosons under the exchange of any pair Xi ^ Xj. In the 
foUowings, we assume our states of Eqs. (1) and (2) have a fixed particle number for 
convenience. It is easy to extend to particle-number nonconserving states as discussed in 
Appendix A. Our final results are valid for both cases because they are expressed in terms 
of density distribution functions. 

Evaluation of kinetic energy involving Jastrow function is helped by Jackson-Feenberg 
transformation [2]. We reproduce the transformation as Eq. (B7) in Appendix B for our 
wavefunctions of Eqs. (33), 



Another equivalent expression is also derived as Eq. (B8), 



(34) 



V? + ^(V^f/) + ^(V.f/)-V, 



(35) 



Both transformations involve one- and two-body density distribution functions only, and the 
Jastrow factor appears on the right of the derivatives. The biased distribution functions 
defined in Eq. (25) are then applicable. The expectation value of a general Hamiltonian 
with an external field (and/or chemical potential) A{x) is calculated as 

fc2 



E 



1 



Using Eq. (34), we derive 



(36) 



El 



dx 



2m 



p[{x) + Aeff{x)pi{x) 



+ - I dXidX2Veff{Xi,X2)p2, 



(37) 



where pi and p2 are the one- and two-body density distribution functions; Aeff{x) is the 
effective external field operator 



Aeff = A{X) - —{V 

8m 



c\2 



(38) 



with operator defined as applying to the Coester states only; Veff is the effective potential 
defined as 



VeffiXi, X2) = v{Xi, X2) - + 'Vl)u{Xi, X2)] 

8m 



(39) 
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and finally p\ is one-body density function derived from the first term of Eq. (34) and written 
in the second quantization form as 



(40) 



Using Eq. (35), we derive second equivalent energy functional as 



E<), 



dx 



—p[{x) + A{x)pi{x) 



+ 



\jdx. 



dxo 



Veff{xi,X2)p2 - ^P2(«) 



(41) 



where p'2{u) is two-body density functions derived from the third integrals of Eq. (35), 



P2(«) 



^'^|e^°r2(«)|*'=). 



with the two-body operator T2{u) given by 



(42) 



T2{u) = lP\xi)lP\x2) [Viu{xi, X2) ■ Vi + V2U{XI, X2) ■ V2] 1p{x2)lp{xi). 



(43) 



The difference between the two energy functionals is that in Ei of Eq. (37) we need to 
take care of the operator which apphes only to the Coester states and in E2 of Eq. (41) 
we need to calculate the two-body density function P2{u). We hope to get experience in real 
applications as which form is more practical. The density distribution functions pi and p2 
were calculated earlier by Eqs. (30) and (32). It is easy to show that p[ and ^2(1*) can be 
calculated in similar fashion. We hence write 

1 



1 



Pi 



(44) 



where fi(x) = e-^T^e^ = 'ij}'' {x)V^'ii}{x) . After evaluating fi(x)|$) = [X[oiF;x) + 
J:iX[j{F;x)C\]\(^) as before, we have 



where 



p[{x)^pf{x) + J2X[4F;x)-^^^ 



pf(x) = {T,r = Xl,(F;x) + j:Xlj{F;x)ri 



(45) 



(46) 



is the similar density function of the Coester states, and the derivative dW^/ dFj in Eq. (45) 
is given by Eq. (27). In similar fashion, we derive 



p'^{u] xi, X2) = P2 + 



(47) 
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where 

p'^{u- xi, X2) = {T^r = X^,o + E Ki9i, (48) 

and X^ Q and Xl^ j are obtained by evaluating f2{u)\^) = (X2 + T,i X2^jCl)\^). 

Eqs. (45) and (47), together with density distributions of Eq. (30) and (32), are all we 
need for calculating the two equivalent energy functional of Eq. (37) and (41). These are 
our main results in this paper. Denoting the three terms from the first term of p[, pi and p2 
as e, and the higher-order remainders as Ai?; with / = 1, 2 for the two energy functionals, 
we rewrite the energy functionals of Eqs. (37) and (41) as 

Ei{F,~g\u) = e{F,~g\u) + m{Frg\u), 1 = 1,2 (49) 

where e(F, g'^, u) is given by, 

e = Kl + J dxA{x)pl{x) + dXidX2Veff{xi,X2)e''^'''''''^pl2ixi,X2), (50) 

with = -{h? /2m){TiY for the kinetic energy of the VCCM states. We notice that 
both energy functions give the same first order approximation e as the term containing 
{y'^Yp\{x) in El vanishes after integration by parts. Furthermore, this first-order energy 
functional e is nothing but the energy functional of the VCCM after replacing the bare 
potential v —>■ V = Veffc'^, namely 



e{F,r,u) = E%F,r)l^y, V{v,u) 



v{xi, X2) - 3— (Vi + Vl)u{xi, X2) 
8m 



(51) 

This is convenient indeed as no new calculations are needed after the VCCM calculations 
have been done. We want to emphasize that, in addition to this simple, intuitively appealing 
approximation of Eq. (51), our main purpose here is to provide the detailed formulas for 
calculating the higher-order terms in AEi of Eq. (49) in the forms of four equations for pi, 
p[, P2 and P2 respectively. 



IV. SUMMARY AND DISCUSSION 



In summary, we have calculated the Hamiltonian expectation value of a general quantum 
many-body system as energy functionals of the Coester-state coefficients {-F/}, bare distribu- 
tion functions {gj} of the Coester states and the Jastrow correlation function u{xi, ^2). Two 
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equivalent energy functionals are derive by Eqs. (37) and (41). A single simple, first-order 
energy from both expressions is derived as the usual VCCM energy functional but with a new 
potential V{v, u) = VeffC^. The formulas for practical calculations of the higher-order terms 
are provided in details. It is easy to see that, due to the short-ranged Jastrow factor, our 
energy functionals do not suffer the divergence as in most coupled-cluster calculations when 
potential w(xi,X2) approaching hardcore potentials. We also want to emphasize that an as- 
sumption has been made in applying the generalized linked- cluster expansion of Eq. (13) for 
the Coester states. We have proved this is valid for a SUB2 truncation in the Coester states 
of spin lattice application (and similar approximations for Bose gas and BCS state), we need 
to examine this validity for the Coester states in a specific truncation approximation for a 
real application. We believe this will not pose major difficulty as we have the simple relations 
between the full distribution functions and the bare distribution functions of Eqs. (23). We 
also like to point out that the generalized linked-cluster expansion of Eq. (13) can not be 
applied to the traditional CCM due to the linear construction of its bra state. 



A similar wavefunctions to 
correlations in nuclear matter 



(33) were employed by Owen for study of spin-dependent 
. In particular, \E'^ was approximated by the product of an 



independent pair function (spin-dependent) and the Slater determinant; and in the cluster 
expansion, the Jastrow function and the independent pair function are treated together. 
This differs from our approach. Our calculations involving Coester states are always in 
the second quantization form and are applicable to higher-order truncation approximations. 
As mentioned in Introduction, KKZ [18| employed the Coester ket state of Eq. (1), and 
their energy and the equations for correlation coefficients were obtained by the Hamiltonian 
eigenequation. The validity of this procedure may be questionable as the Coester ket state 
was approximated in a truncation approximation. In the traditional CCM approach, the 
use of Hamiltonian eigenequation is completely equivalent in any truncation approximation 
;o a bi-orthogonal variational equations where the specific bra state of Eq. (7) is employed 
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24| . due to the fact that the normalization integral is always unity and that the bra 



state is linear in the coefficients. This is certainly not the case after inclusion of the Jastrow 
function. In particular, the normalization integral in KKZ's equations must contain an 
additional term involving the bra-state coefficients of Eq. (7). It is not clear how this change 
will complicate their analysis. However, we agree with KKZ and believe that the Jastrow 
theory and coupled-cluster theory in principle complement each other and combination of 
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these two theories may overcome each other's difficulties. The approach presented here is our 
such attempt and we have taken the advantage of the recent progress in the coupled-cluster 
theory. 

It is interesting to compare our energy functional E{F,g'^, u) with the counterpart in the 
traditional Jastrow theory, E{a,u), where a is the one-body function. Clearly, the missing 
state-dependent correlations in E{a,u) are now included in E{F,g'^,u) in terms of F and 
g'^. However, a typical calculation of the Jastrow theory practiced today mostly includes 
resummation of all cluster terms of the linked cluster expansion by the hyper-netted chain 
(hnc) approximation for bosons or Fermi-hnc approximation for fermions j^, and one 
optimization route for the boson system is provided by the pair-phonon analysis (PPA) of 
Campbell and Feenberg [32'|. This is possible as the reference states are single-particle states. 
In our cluster expansion calculations, the reference state is the Coester states which already 
contain rich correlations including in particular correct long-ranged correlations such as in 
the SUB2 approximation. The introduction of the Jastrow function is to provide correct 
description of strong, short-ranged correlations and we expect such scheme may provide 
reasonable results even if we include a ffist few cluster terms. It will also be interesting to 
investigate the relations between our approach and the PPA of Campbell and Feenberg. 
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APPENDIX A: A-REPRESENTATION AND GENERALIZED LINKED- 
CLUSTER EXPANSION 



We first consider the conversion relations between real-space wavefunctions and the cor- 
responding states in momentum-space. Let \xi, ■ ■ • ,xn) be the single-particle basis of a 
A^-particle system, namely 

• • • ,XAr) = ^V'^M • • • V'^(xi)IO) (Al) 

where ip'^{x) are the field operators and |0) is the corresponding vacuum state. The unity 
operator is then given by 

l^Z] / dX\xi, - ■ ■ ,xn){xi, - ■ ■ ,xn\, (A2) 

where / dX — J dxi ■ ■ ■ dxjq. The inner product between two arbitrary states can be written 
as 

(^il^r^) = ^ I dX(*i|xi,---,x^)(xi,---,x^|*2) = E / (A3) 

where = • • • , xjv) = (xi, • • • , xa^I^'i) and ^'2 = ^'2(0:1, • • • , xjv) = • • • , xjv|*2) 

are the corresponding wavefunctions in real space with proper boson or fermion symmetry. 
If these two states have a fixed particle number N, the above inner product is then simply 
given by the familiar quantum mechanics formula 

(*i|*2) = / (i^*t*2, (A4) 

because {xi, ■ ■ ■ ,xiyi\'^2) = ^n',n'^2{xi, ■ ■ ■ ,Xn). In the foUowings, wc assume our states 
have a fixed particle number for convenience. It is easy to extend to particle-number 
nonconserving states by using Eq. (A3). Our final results in this paper are vahd for both 
cases because they are expressed in terms of density distribution functions. 

We denote ^'^(xi, • • • , xn) = (xi, • • • , xnI"^"^) as the Coester state function in real space, 
where Coester state l^'^) is given by Eq. (4). The application of the Jastrow operator S*" of 
Eq. (3) on the single-particle basis can be written as 
go u 

— \xi,---,xn) = —\xi,---,xn), U = ^u{xi,Xj) (A5) 

or 

e^°/Vi,---,a;jv) = e^/Vi,---,a;iv). (A6) 
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We therefore have, usina; the fact that real function, 

vl/(xi, • • • , x^) = (xi, • • • , XN\e'°/^m = e^'/'ix,, • • • , x^vl^^) = ^"^^ (A7) 

where = e^/^ = exp {j2i<j u{xi, Xj) is the famihar Jastrow wavefunction. Similarly, 
the bra-state wavefunction is written as 

^(xi, ■■■,xn) = = ^"e^/^ (A8) 

where are the corresponding bra-state wavefunctions of Eqs. (6) or (7) in real space. The 
Coester wavefunctions \E''^ and obey proper symmetry, namely they are antisymmetric 
for fermions and symmetric for bosons under the exchange of any pair Xi ^ Xj. We do not 
need to know their explicit functional forms as our later calculations are always carried out 
in second quantization form of momentum space. 

We follow the similar analysis for the evaluation of Eg. (13) as in the traditional Jastrow 
theory [2, Il6|]. After the usual cluster expansion of the Jastrow wavefunction in terms of 
Yvon-Mayor functions Zn{Y) with the bound function Y12 = e"^^^'^^^ — 1, 

= 1 + ^Z2{xi,Xj) + Yj Z3{xi,Xj,Xk) ^ , (A9) 

i<j i<j<k 

the expectation in Eq. (13) is written as, 

{eS°y = 1+ J2 7T7 '-\n dXl--- dXN^'ZrriXi, ■ ■ ■ , Xn)-^" = 1 + ^ ^ dXi ■ ■ ■ dXnZnP^ 

(AlO) 

where Ic = is the normalization integral, the first few Zn are given by 

Z2 = yi2, Z3 = Y12Y23 + 123^31 + YsiYu + ^12^23^31; (All) 

etc., and is the n-body density distribution functions 

Pn = ■ ■ ■ V^t(x„)^(x„) ■ ■ ■ ^{x^)y. (A12) 

Evaluation of of the Coester states in terms of bare distribution function g'^ were discussed 
in details by Eqs. (18-23) of Sec. II. Here we need to consider their cluster properties in 
order to apply the linked-cluster theorem. For this purpose, we introduce the so-called A- 
representation discussed as follows. We consider evaluation of n-body distribution functions 
in the state with the following cluster decomposition property: any higher-order can 
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be written as products of lower order with m — 1,2, • ■ ■ ,n — 1, plus a nondecomposable 
core A„. Let us demonstrate this property in details. We start with the one-body density 
distribution matrix, also the one-body core by definition, 

pf{x; x') = {i^\x)i^{x')Y = Ai(x; x'), (A13) 

where we have used the usual notation (■•■)"' for the expectation value in state 1^*^). The 
two-body distribution (diagonal) function, p2 — {xi)ip^ {x2)ip{x2)ip{xi)y is then given by 
using this decomposition property, 

P2 = Ai(a;i; a;i)Ai(x2; X2) - Ai{xi; a;2)Ai(a;2; xi) + A2(a;i, 0:2; xi, X2), (A14) 

where the first two terms are as given before for the fermions (for bosons, all terms have 
positive sign) and the A2 is the 2-body nondecomposable core. We can also write Eq. (A14) 
in a symbolic notation as 

= Ai * Ai A2. (A15) 
In similar way, the three-body density distribution function is then given by 

= Ai * Ai * Ai Ai * A2 A3, (A16) 

where A3 is the nondecomposable core. All terms from the product Ai * Ai* Ai of Eq. (A16) 
were given in details as diagrams in the references quoted earlier. Similar to these dia- 
grams, the new contributions in products Ai * A2 also include both unlinked terms such as 
Ai(xi; xi)Ai(x2, X3; X2, X3), and linked terms such as Ai(xi; X3)A2(x2, X3; X2, xi). The rules 
for sign and for the symbohc product (*) can also all be defined in similar fashion as before 
and wc will discuss them in details somewhere else. Such decomposition can be carried to 
higher-order. 

Using the cluster decomposition property as discussed above, we have the key ingredient 
for the linked-cluster expansion, namely, all contributions in the product Znp^ can be repre- 
sented by clusters of diagrams. If we denote linked diagrams as Fa, Fb, etc., a contribution 
in similar cluster expansion as Eq. (AlO) for state |\1''^) can then be written as 

r = (r^)^^(rB)^----, a^b, (ai7) 

with coefficients 

N{A,B,---) = ^ ,^ • ,^ (A18) 
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giving by the number of distributing the hnked part along the n points of the diagram. We 
therefore have 

EV"= E -,N{A,B,---){TAr^{TBr--- = eM'^A/nA\ + TB/nB\ + ■■■). (A19) 

Hence we have the following generalized linked-cluster expansion for the state l^l/'^), 

ln(e^")" = ^ 1 /■ dxi ■ ■ ■ rfx„[Z„p^]i, (A20) 

n=2 •' 

where the notation [Znpnjh denotes the contributions to ZnPn by linked diagrams only. Since 
all terms in Z2 and are linked, the first two terms in Eq. (A20) are independent of the 
diagram structures of p2 and pf and we write 

ln(e^°)°' = ^ / dxidx2Z2pi + ^ / dxidx2dx3Z3p^ + E ~T f dxi ■ ■ ■ dxn[Znpt]L- (A21) 

A J J n=4 

From n = 4 and onward, we need to know the diagram details of for calculating their 
contributions. 

In order to apply the generalized linked-cluster expansion of Eq. (A21) to our Coester 
states, we need to prove that the Coester states satisfy the cluster decomposition property 
as discussed above. In our earlier VCCM calculation for spin lattices, we have shown in- 
deed the Coester states satisfy such property in a SUB2 approximation employed, where 
arbitrary order distribution functions can be calculated by the simple functional derivative 
Qiiji /dFij = gijiQiij and these bare distribution functions correspond to the density distribu- 
tion matrices (similar analysis also applied to the SUB2 state for the Bose gas and the BCS 



superconductors) [20|, l27|] . We also notice that a similar so-called SUBm truncation approx- 
imation in the A-representation can also be defined as the approximation retaining up to m 
core distribution tensors only. It is intuitive to relate the real-space cluster parametrization 
by core distribution tensor {A„} in the A-representation and momentum-space parametriza- 
tion by {Fi.gi} in the Coester states; the Coester representation provides a practical way 
to calculate these core tensors. We will not intend to provide a general proof that the Co- 
ester states in any truncation approximation will satisfy the cluster decomposition property. 
We will adopt a practical strategy and apply the the linked-cluster expansion formula of 
Eq. (A21) to the Coester states in real applications and examine the cluster property in the 
particular truncation approximation employed. We believe this will not cause a major diffi- 
culty as we use the relation between full distribution functions and the bare distribution 
functions gj as given by Eqs. (23). 
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We also like to point out that the traditional CCM states certainly fail the cluster de- 
composition property due to the linear construction of the bra state of Eq. (7) . This can be 
easily seen as any expectation in the CCM is always linear in the bra state coefficients, con- 
tradictory to the cluster decomposition property. Therefore, the generalized linked- cluster 
expansion of Eq. (A21) can not be applied to the CCM states. 



APPENDIX B: JACKSON-FEENBERG TRANSFORMATION 



We next consider the expectation value of the kinetic energy operator using wavefunctions 
of Eqs. (A7) and (A8). We follow the derivation as given in Ref. 2 but keep using the notation 
as it is different to the ket-state counterpart in the traditional CCM. Applying the nested 
commutation formula in the following integral 



we have 



j dX^'^e^ 

We can also apply the same nested commutation formula differently as. 



(Bl) 



(B2) 



and derive another equivalent expression as 

1 



dX<i/' 



V,:--(V.t/) 



(B3) 



Summing Eqs. (B2) and (B3), the two (VjC/)^ terms (3-body terms) cancel each other out, 
we have 

+ [e^(V,C/)V, + V,e^(V,f/)] ^"j, 



and integrating by parts for the third and last terms, we have 



2L 



-I 



dX 



where we have introduced Bi as 

1 



S,: = - / dX 



(B4) 



(B5) 
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If — these two terms cancel each other out and Bi — 0, we therefore have 



dX 



Using the identity 



the well-known Jackson-Feenberg formula is then derived as, with = ^ 



(B6) 



(B7) 



(B8) 



Eq. (B6) can also be written as after integrating by part for the last term, 
li^ J dX^%^ 

We notice that in Eqs. (B7) and (B8) we have assumed the Coester ket- and bra-state 
wavefunctions equal, = \E''^. This is true for the VCCM states of Eq. (4-6) but not true 
in general for the CCM state of Eq. (7) in any finite truncation in 5" and 5". Therefore the 
two terms in of Eq. (B5) do not cancel in general in the the CCM basis. Our conclusion 
is that the kinetic expectation value in the CCM basis will also contain three-body density 
distribution functions. 
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